7 induced multiparticle emissions of medium mass nuclei at intermediate energies 
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A comprehensive analysis of multiparticle emissions following photon induced reactions at inter- 
mediate energies is provided. Photon induced reaction is described in the energy range of ~ 30 — 140 
MeV with an approach based on the quasideuteron nuclear photoabsorption model followed by the 
process of competition between light particle evaporation and fission for the excited nucleus. The 
evaporation-fission process of the compound nucleus is simulated in a Monte-Carlo framework. The 
study shows almost no fission events for the medium mass nuclei and reproduces the available 
experimental data of photonuclear reaction cross sections satisfactorily at energies ~ 30 — 140 MeV. 
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With energetic incident photon it is possible to induce 
nuclear reactions including fission in most elements. Pho- 
tonuclear reactions in different energy ranges, like giant 
dipole resonance (GDR) and quasideuteron (QD) energy 
regions, have been studied in the past as they provide a 
wide range of information either on the initial nuclear ex- 
citation mechanism or on characteristics of the compound 
nucleus decay channels. In fact photoneutron cross sec- 
tions in the GDR energy region have been compiled [l[ for 
most of the nuclei in the periodic table, while measure- 
ments in the QD region were mainly focussed on heavy 
nuclei, thus leading to a lack of study on intermediate 
nuclei. For heavier nuclei, particularly for actinides and 
preactinides, with high enough energy of incident pho- 
tons, the dominant reaction mechanism is fission where 
photofission cross sections are quite large 0]. With de- 
creasing fissility (Z 2 /A) of the target nucleus and energy 
of incident photon, fission probability decreases, whereas 
other reaction channels such as (7, xn) and (7, xnyp) type 
become more important. Thus it becomes important to 
study the photonuclear reactions in the intermediate en- 
ergy range for medium mass elements. Investigations 
of such multiparticle emissions following photon induced 
reactions in the intermediate energy region in medium 
mass elements may reveal important features of the na- 
ture of the photoabsorption mechanism as well as the 
decay channel characteristics. 

The aim of the present work is to investigate photonu- 
clear reactions of the type (7, xn.) and (-y,xnyp). The 
present study is restricted to the QD region which is 
above the GDR or the isovector giant quadrupole reso- 
nance (IVGQR) regions and is below the pion threshold. 
The different photonuclear reactions for V, Zn, Sn, Ce, 
Sm, Yb, Ta, Au, Pb and Bi are consistently described as 
a two step process. In the rapid stage of a photonuclear 
reaction, the incoming 7 is assumed to be absorbed by a 
neutron-proton [n-p] pair inside the nucleus. The rapid 



stage is then followed by a subsequent de-excitation of the 
compound nucleus (CN) via evaporation. Such a statis- 
tical decay of the CN is the slow stage of a photonuclear 
reaction. The quantitative description of the process is 
based on the liquid drop model (LDM) for nuclear fission 
by Bohr and Wheeler [3[ and the statistical model of nu- 
clear evaporation developed by Weisskopf 4]. For this 
slow stage, an evaporation Monte-Carlo routine, based 
on the Weisskopf statistical theory, is used to address 
the de-excitation of the CN in terms of the competition 
between particle evaporation and nuclear fission. 

The dominant mechanism for nuclear photoabsorption 
at intermediate energies is described by the QD model 
H which is employed to access [6| the total photoabsorp- 
tion cross section in nuclei. It is based on the assumption 
that the incident photon is absorbed by a correlated n-p 
pair inside the nucleus, leaving the remaining nucleons as 
spectators. Such an assumption is enforced when wave- 
length of the incident photon is small compared to nu- 
clear dimensions but not too small, where the wavelength 
becomes smaller than the average internucleon spacing in 
the nucleus. The total nuclear photoabsorption cross sec- 
tion erj is then proportional to the available number of 
n-p pairs inside the nucleus and also to the free deuteron 
photodisintegration cross section ad(E 7 ) : 



(1) 



where N, Z and A are the neutron, proton and mass 
numbers respectively, 4 factor represents the fraction of 
correlated n-p pairs and the function e~ D i E ~i accounts 
for the reduction of the n-p phase space due to the Pauli 
exclusion principle. A systematic study of total nuclear 
photoabsorption cross section data in the intermediate 
energy range shows that D = 0.72A 081 MeVQ. The free 
deuteron photodisintegration cross section |8] is given by 
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where E* = E* — B n and E* f = E* — Bf are the nu- 
clear excitation energies after the emission of a neutron 
and after fission, respectively, where B n is the binding 
energy of the emitted neutron. T n and T f are the partial 
widths for the decay of the excited compound nucleus 
via neutron emission and fission, respectively, and the 
parameters a n and af are the level density parameters 
for the neutron emission and the fission, respectively and 
A'o = h 2 /2mrQ where m and ro are the neutron mass and 
radius parameter respectively. The emission probability 
of particle k relative to neutron emission is calculated 
according to the Weisskopfs statistical model [J] 
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where B (= 2.224 MeV) is the binding energy of the 
deuteron. The QD model H [f| of nuclear photoabsorp- 
tion is used together with modern rms radius data to ob- 
tain Levinger's constant L — 6.8— 11.2A~ 2 / 3 + 5.7 A~ 4 / 3 
of nuclei throughout the Periodic Table and is in good 
agreement Q with those obtained from the experimen- 
tally measured values. 

At the QD energy range, as a consequence of the 2? 
primary photointeraction, 7+(n+p) — > n*+p*, in most E, 
of the cases excited compound nucleus is formed with 
the same composition as the target nucleus where both 
neutron and proton are retained inside the nucleus and 
the probabilities that either the neutron or the pro- 
ton or both escape from the nucleus are extremely low. 
Hence, the recoiling nucleus can be viewed as a com- 
pound nucleus having the same composition as the tar- 
get nucleus but with excitation energy E* = moC 2 [(l + 
2A 7 /m c 2 ) 1 / 2 — 1] where E 1 is the incident photon 
enegy and m is the rest mass of the nucleus before pho- 
ton absorption. This excited compound nucleus then un- 
dergoes successive evaporation of neutrons, protons and 
light particles or fission. Hence the photonuclear reaction 
cross section a r is a product of the nuclear photoabsorp- 
tion cross section <rj and the statistical decay probabil- 
ity p for a reaction channel and is, therefore, given by 
a r (E 7 ) — a^(E 7 )p. The decay probability p is basically 
the branching ratio T r /T where T r and T are the partial 
and total reaction widths respectively. 

The evaporation stage is calculated by applying the 
statistical theory 0] . The decay of the CN takes into ac- 
count all accessible channels, with the related branchings 
calculated in terms of the nuclear level densities of the 
daughter nuclei. The basic steps are the calculations of 
the relative probabilities between the competing channels 
(particle evaporation and nuclear fission). The probabil- 
ity of fission relative to neutron emission is calculated 
using Vandenbosch-Huizenga's equation [ll|, given by 
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FIG. 1: The plots of cross sections oer as a function of mass 
number An of the evaporation residues resulting from (7, xn) 
and (7, xnyp) types of reactions for 64 Zn at E 1 — 30, 40 MeV. 
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where is the level density parameter for the emis- 
sion of the particle k, Jk/jn = 1 for k — p, 2 for 
k = 4 Hc, I for k = 2 H, 3 for k = 3 H and 2 for k = 3 Hc. 
El = E* — (Bf. + Vfe) is the nuclear excitation energy 
after the emission of particle k [T3] • Bk are the binding 
energies of the emitted particles and are the Coulomb 
potentials. The evaporation- fission competition of the 
compound nucleus is then described in the framework 
of a projection angular momentum coupled [l3j Monte- 
Carlo routine [LH ]. Any particular reaction channel r is 
then defined as the formation of the compound nucleus 
via photoabsorption and its decay via particle emission 
or fission. Thus, fission is considered as a decay mode. 
The photonuclear reaction cross sections a r are calcu- 
lated using the equation a r — a^n r /N where n r is the 
number of events in a particular reaction channel r and 
N is total number of events that is the number of the 
incident photons. 

Each calculation is performed with 40000 events us- 
ing a Monte-Carlo technique for the evaporation-fission 
calculation. This provides a reasonably good computa- 
tional statistics. However, possible contributions from 
the direct reactions are omitted from the present cal- 
culations. The cross section a r (E 1 ) for any particular 
reaction channel r is then calculated using the equation 
a r (E 1 ) = erj \E~ ( )n r / 'N . Detailed calculations have been 
performed for 51 V, 64 Zn, 118 Sn, 140 Ce, 154 Sm, 174 Yb, 
181 Ta, 197 Au, 208 Pb and 209 Bi at E 1 = 30 MeV to 140 
MeV. In Table 1 the results for {j,xn) and (7, xnyp) re- 
action cross sections are provided. The statistical error 
in the theoretical calculations for the photonuclear re- 
action cross sections can be estimated using the equa- 
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TABLE I: Photofission cross section and three largest photonuclear reaction cross sections in (7, xnyp) decay channels. The 
number of evaporated neutrons x and protons y are provided as (x, y) adjacent to each reaction cross sections. 



Target Z 2 /A 


Calculated 


E-y [MeV] 


£ 7 [MeV] 


E-y [MeV] 


E~ f [MeV] 


E~ f [MeV] 


E-y [MeV] 


E-, [MeV] 


nuclei 




quantity [mb] 


30 


40 


60 


80 


100 


120 


140 


51 V 

V 


1 n V7 


O a 


14.11 


10.92 


7.07 


5.01 


3.78 


2.97 


2.42 






a r (-y,xnyp) 


10.76 (2,0) 


4.30 (3,0) 


3.46 (3,1) 


1.80 (4,1) 


0.95 (4,2) 


0.55 (5,2) 


0.31 (5,3) 








2.72 (1,1) 


2.36 (2,0) 


1.12 (4,0) 


0.87 (3,2) 


0.58 (5,2) 


0.38 (5,3) 


0.25 (7,4) 








0.32 (1,0) 


2.28 (2,1) 


0.52 (2,2) 


0.46 (4,2) 


0.49 (5,1) 


0.26 (4,3) 


0.23 (8,4) 






Of 


0.0 


0.0 


0.18xl0" 3 


0.25xl0" 3 


0.85xl0 -3 


1.71xl0" 3 


2.48xl0" 3 


64 Zn 


14.06 


„T 
O a 


16.13 


12.84 


8.57 


6.16 


4.68 


3.71 


3.03 






a r (-y,xnyp) 


9.67 (1,1) 


4.67 (2,1) 


3.19 (2,2) 


1.64 (3,2) 


0.82 (4,3) 


0.54 (6,5) 


0.43 (5,4) 








2.73 (3,2) 


2.78 (1,2) 


1.83 (4,3) 


1.00 (5,3) 


0.79 (5,4) 


0.52 (5,3) 


0.39 (6,5) 








2.08 (2,0) 


1.42 (3,2) 


1.03 (3,1) 


0.67 (4,4) 


0.50 (4,2) 


0.40 (6,4) 


0.32 (6,6) 






a f 


0.0 


0.0 


0.0 


0.15xl0 -3 


0.0 


0.28xl0 -3 


0.45xl0 -3 


118 Sn 


21.19 


o a 


19.32 


17.21 


12.84 


9.76 


7.67 


6.22 


5.16 






a r (-y,xnyp) 


12.95 (2,0) 


15.03 (3,0) 


9.44 (5,0) 


5.74 (6,0) 


2.93 (8,0) 


1.63 (9,0) 


1.28 (9,1) 








6.33 (3,0) 


1.91 (4,0) 


2.57 (4,0) 


2.12 (7,0) 


2.03 (7,0) 


1.57 (8,1) 


0.97 (10,0) 








0.04 (1,1) 


0.14 (2,1) 


0.34 (4,1) 


0.86 (5,1) 


1.06 (7,1) 


0.69 (8,0) 


0.69 (10,0) 






a f 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


140 Ce 


24.03 


o a 


19.37 


18.00 


14.01 


10.88 


8.66 


7.08 


5.91 






a r (-y,xnyp) 


16.06 (2,0) 


15.96 (3,0) 11.35 (5,0) 


4.81 (6,0) 


4.19 (8,0) 


2.14 (9,0) 


1.44 (10,1) 








3.23 (3,0) 


1.66 (4,0) 


1.48 (4,0) 


3.94 (7,0) 


1.45 (7,1) 


1.43 (8,1) 


1.03 (9,1) 








0.06 (1,1) 


0.25 (2,1) 


0.58 (4,1) 


1.04 (5,1) 


1.15 (7,0) 


0.94 (10,0) 0.72 (10,0) 






Of 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


154 Sm 


24.96 


o a 


19.07 


18.20 


14.54 


11.44 


9.19 


7.55 


6.32 






°r('j,xnyp) 


17.82 (3,0) 


17.05 (4,0) 11.44 (6,0) 


7.62 (8,0) 


4.66 (10,0) 3.54 (11,0) 2.50 (12,0) 








0.83 (2,0) 


0.84 (3,0) 


1.61 (7,0) 


2.04 (7,0) 


3.05 (9,0) 


1.45 (10,0) 0.84 (13,0) 








0.42 (4,0) 


0.28 (5,0) 


1.28 (5,0) 


1.12 (9,0) 


0.42 (10,2) 0.76 (12,0) 0.64 (11,1) 






Of 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


174y b 


28.16 


„T 
O a 


18.67 


18.49 


15.33 


12.29 


9.97 


8.25 


6.95 






a r (-y,xnyp) 


18.02 (3,0) 


17.38 (4,0) 11.72 (6,0) 


5.79 (8,0) 


5.65 (9,0) 


3.27 (11,0) 2.92 (12,0) 








0.62 (2,0) 


1.04 (3,0) 


3.34 (5,0) 


5.64 (7,0) 


1.65 (8,0) 


2.95 (10,0) 0.98 (11,0) 








0.03 (4,0) 


0.05 (5,0) 


0.08 (4,1) 


0.30 (6,0) 


1.64 (10,0) 


0.44 (9,1) 


0.91(13,0) 






Of 


0.0 


0.0 


0.0 


0.0 


0.0 


0.62xl0" 3 


0.0 


181 Ta 


29.44 


„T 
O a 


18.50 


18.55 


15.58 


12.57 


10.24 


8.49 


7.17 






a r (-y,xnyp) 


17.57 (3,0) 


17.18 (4,0) 11.68 (6,0) 


5.92 (7,0) 


6.16 (9,0) 


3.24 (11,0) 2.80 (12,0) 








0.92 (2,0) 


1.32 (3,0) 


3.56 (5,0) 


5.68 (8,0) 


1.39 (10,0) 2.77 (10,0) 0.91 (11,1) 








0.01 (1,1) 


0.01 (4,2) 


0.08 (4,1) 


0.29 (6,1) 


1.34 (8,0) 


0.59 (9,1) 


0.89 (11,0) 






Of 


0.0 


0.0 


0.0 


3.1xl0" 3 


1.02xl0" 3 


1.7xl0" 3 


3.1xl0" 3 


197 Au 


31.68 


„T 
O a 


17.97 


18.55 


16.03 


13.12 


10.78 


8.99 


7.62 






a r (^,xnyp) 


15.70 (3,0) 


14.04 (4,0) 10.10 (6,0) 


7.28 (7,0) 


6.44 (9,0) 


3.89 (10,0) 3.87 (12,0) 








2.28 (2,0) 


4.49 (3,0) 


5.77 (5,0) 


4.91 (8,0) 


2.28 (8,0) 


3.27 (11,0) 0.86 (11,0) 








0.9xl0~ 3 (1,1) 


0.01 (2,1) 


0.06 (4,1) 


0.57 (6,0) 


1.25 (10,0) 0.40 (10,1) 0.80 (13,0) 






Of 


0.0 


1.85xl0" 3 


9.22xl0" 3 


0.02 


0.05 


0.07 


0.16 


208 Pb 


32.33 


o a 


17.48 


18.39 


16.21 


13.39 


11.07 


9.27 


7.88 






G r {^,xnyp) 


15.47 (3,0) 


16.12 (4,0) 11.79 (6,0) 


6.49 (8,0) 


6.97 (9,0) 


5.12 (11,0) 3.53 (12,0) 








2.01 (2,0) 


2.23 (3,0) 


4.27 (5,0) 


6.46 (7,0) 


2.41 (10,0) 2.61 (10,0) 2.16 (13,0) 








5.68 xl0~ 3 (4,0) 


0.01 (5,0) 


0.03 (4,1) 


0.19 (6,0) 


1.18 (8,0) 


0.53 (12,0) 0.77 (11,0) 






Of 


1.75 xl0~ 3 


0.02 


0.06 


0.12 


0.20 


0.33 


0.44 


209 B j 


32.96 


Oa 


17.49 


18.43 


16.28 


13.46 


11.13 


9.33 


7.93 






a r (-y,xnyp) 


15.71 (3,0) 


15.91 (4,0) 10.90 (6,0) 


6.65 (7,0) 


6.60 (9,0) 


3.97 (11,0) 2.40 (12,0) 








1.75 (2,0) 


2.42 (3,0) 


4.86 (5,0) 


5.57 (8,0) 


1.72 (10,0) 2.57 (10,0) 1.36 (13,0) 








6.12xl0~ 3 (2,1) 


0.03 (2,1) 


0.14 (4,1) 


0.33 (6,1) 


1.11 (8,0) 


0.57 (10,1) 0.93 (11,1) 






Of 


0.02 


0.06 


0.21 


0.40 


0.60 


0.88 


1.18 
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FIG. 2: Plots of total photoabsorption cross section versus 
photon energy. The continuous lines represent the total pho- 
toabsorption cross sections from the present calculation for 
discrete isotopes 118 Sn, 140 Ce, 181 Ta and 208 Pb whereas the 
data points are those corresponding to the measured values 
[T5I ] for natural targets of Sn, Ce, Ta and Pb. 
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FIG. 3: Plots of average neutron multiplicities versus photon 
energy. The continuous lines represent the average neutron 
multiplicities from the present calculation for discrete isotopes 
118 Sn, 140 Ce, 181 Ta and 208 Pb whereas the data points are 
those corresponding to the measured values [l6j] for natural 
targets of Sn, Ce, Ta and Pb. 



tion a r ± Aoy = cj [?v ± ^/n^]/N which implies that 
Aoy = \J aJa r /N. In Fig.l the plots of cross sec- 
tions of evaporation residues resulting from (7, n) and 
(7, np) type of reactions for 64 Zn for E 1 = 30 MeV 
and 40 MeV are shown as typical cases. In Fig. 2 the 
plots of total photoabsorption cross section versus pho- 
ton energy E 7 are presented. The continuous lines rep- 
resent the total photoabsorption cross sections for Sn, 
Ce, Ta and Pb from the present calculations whereas 
the points with error bars are the experimental data 
[l5T| , Fig. 3 shows the plots of average neutron multiplici- 
ties < v >=T, x xa('j, xnyp)/<Ja versus photon energy E 1 . 
The continuous lines represent the average neutron mul- 
tiplicities for Sn, Ce, Ta and Pb from the present calcu- 
lations whereas the points with error bars are the exper- 
imental data [TB|. Considering the fact that the present 
calculations for photoabsorption cross sections are for 
discrete isotopes 118 Sn, 140 Ce, 181 Ta and 208 Pb whereas 
the measurements fl5| were performed using natural tar- 
gets of Sn, Ce, Ta and Pb, present results agree quite 
closely with the available experimental data at interme- 
diate energies. This effect is more pronounced in the plots 
for < v > where at E 1 above 50 MeV, present calcula- 
tions gradually overestimate. This behaviour is expected 
as neutron yields would be less in presence of lower mass 
isotopes in natural targets. Although the present simple 
model appears to start breaking down above photon en- 
ergies 50-60 MeV, more experimental data are required 
before any definitive conclusion can be reached. It is 
worthwhile to mention here that the present formalism 
provides excellent estimates [2] for photofission cross sec- 
tions of actinide and pre-actinide nuclei at intermediate 
energies. 

The agreement of the present theoretical calculations 
according to the QD model with the experimental data 
for photoabsorption cross sections a^(E 7 ) is much bet- 
ter than that for the average neutron multiplicities < v > 
which can be visualized from Fig. 2 and Fig. 3 respectively. 
Therefore, the main physics issue remains in the second 
stage of the process. In an earlier work, good agreement 
for average neutron multiplicities < v > were achieved 
by inclusion of the precompound decay process [l7| to 
account for higher neutron yield than estimated by the 
evaporation model alone. Certainly, similar modification 
can only worsen the present results. The implication of 
such a state of affairs suggests, perhaps, the need for 
more precise measurements since even one or two miss- 
ing neutrons in the measurements can account for the 
discrepancy in < v > . This point is further strengthened 
by a recent work [l8[ where statistical decay of the ex- 
cited compound nucleus was simulated in a Monte-Carlo 
framework like the present work 0] to describe the evap- 
oration process. Both provided excellent estimates for 
photofission cross sections of actinide nuclei in a wide 
range of energies {E 1 = 68 MeV - 3.77 GeV), ranging 
from the QD region (~ 30 — 140 MeV) @ to the region 
of intranuclear cascade model (above the pion production 
threshold) [H|. 
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In summary, the cross sections for the fission and 
the evaporation residues are calculated for photon in- 
duced nuclear reactions at intermediate energies. Monte- 
Carlo calculations for the evaporation-fission competi- 
tion are performed assuming 40000 incident photons for 
each calculation which provides a reasonably good statis- 
tics for computationally stable results. Present calcu- 
lations provide good estimates of cross sections for the 



reaction channels (7, xn) and ( , y,xnyp) for nuclei 51 V, 
64 Zn, 118 Sn, 140 Cc, 154 Sm, 174 Yb, 181 Ta, 197 Au, 208 Pb 
and 209 Bi at E 7 = 30 MeV to 140 MeV. No fission event 
has been observed below 40 MeV for these nuclei except 
few events in cases of 209 Bi and 208 Pb. Additional exper- 
iments may be necessary for more precise measurements 
of neutron multiplicities. 
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